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ABSTRACT 

Aims. The radiative counterpart of the supermassive black hole at the Galactic center (GC), Sgr A*, is subject to frequent flares visible 
simultaneously in X-rays and near-infrared (NIR). Often, enhanced radio variability from centimeter to sub-millimeter wavelengths 
is observed to follow these X-ray/NIR eruptions. We present here a multi-wavelength campaign carried out in April 2009, with the 
aim of characterizing this broadband flaring activity. 

Methods. ConcuiTent data from the XMM-Newton/EPlC (2-10 keV), VLT/NACO (2.1 yum, 3.8 yum), APEX/LABOCA (870 yum), 
and FermifLAT (0.1-200 GeV) instruments are employed to derive light curves and spectral energy distributions of new flares from 
Sgr A*. 

Results. We detected two relatively bright NIR flares both associated with weak X-ray activity, one of which was followed by a strong 
sub-mm outburst ~200 min later. Photometric spectral information on a NIR flare was obtained for the first time with NACO giving a 
power-law photon index a = -0.4±0.3 (Fy oc v"). The first attempt to detect flaring activity from the Fermi GC source IFGL J1745.6- 
2900 is also reported. NIR, X-ray, and sub-mm flares are finally modeled in the context of non-thermal emission processes. It is found 
that the simplest scenario involving a single expanding plasmoid releasing synchrotron NIR/sub-mm and synchrotron self-Compton 
X-ray radiation is inadequate to reproduce the data, but suggestions to reconcile the basic elements of the theory and the observations 
are proposed. 

Key words. Galaxy: center - Black hole physics - Radiation mechanisms: non-thermal - X-rays: general - Infrared: general - 
sub-millimeter - Gamma-rays: observations 



1. Introduction 

The tracing of stellar orbits in the heart of the Milky Way since 
the early 90's has definitely established the presence of a central 
high concentration of mass, ~4 x 10^ Mq, confined to a region 
no bigger than the Solar System ( |Ghez et al.| 2008[ |Gillesseii| 
|et al.| 2009 [ )■ This is the best evidence to date for the existence 
of a supermassive black hole at the Galactic center (GC, see 
|Melia|[2007[ for a comprehensive review). The relative prox- 
imity of this black hole (8 kpc, Reid et al. 2009) has enabled 
thorough observations of its weak point-like radiative manisfes- 
tation, Sagittarius A* (Sgr A*), all accross the electromagnetic 
spectrum. In spite of a high interstellar extinction, wich pre- 
cludes detections at optical and ultraviolet wavelengths, Sgr A* 
has now been positively detected in radio (from cm to sub-mm 
wavelengths, |Balick & Brown||1974[ [ZyUca & Mezger||1988p, 
near-infrared (NIR, Genzel et al.||2003 i, and X-rays (Baganoff 
|et al.|20()3] l. In the vast y-ray domain, point-like sources coinci- 
dent with the position of Sgr A* have been found at the 10 keV 



(IGR J 17456-290 l,|Belanger et al.|200 4i GeV (1FG LJ1745.6- 
290 0,|Abdo et al.|2C)10| l, and TeV (HESS J 1 745-290, |Aharonian| 
et al. 2004Jl levels, though their association with the supermas- 
sive black hole is less evident. 

One striking feature of Sgr A* is its temporal variabil- 
ity. Albeit much more quiet than typical active galactic nuclei 
(AGNs), Sgr A* displays rapid flares from time to time, in X- 
rays and NIR, thus temporarily enhancing its luminosity by 
factors up to ~160 on timescales of 1-2 hr X-ray flares have 
been discovered at a rate of about one event per day by the 
Chandra (Baganoff et al. 2001 ) and XMM-Newton (Goldwurm] 
|et al.|2003) satellites, while NIR flares seem to recur more fre 



quently (2-6 events per day) as observed by the VLT ( Genzel 
et al."2003'), the Keck ( |Ghez et al.|2"004ll, the HST ( |Yusef-Zadeh" 



et al. 2006aj ), and Subaru ( Nishiyama et al.|2009 
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These flares usually present short timescale variations point- 
ing to an emitting region with a size on the order of IQR^ or 
less, where is the Schwarzschild radius of the black hole. 
They have thus offered new probes into the central engine and 
received considerable attention via extensive multi-wavelength 
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campaigns since their discoveries. On the one hand, such cam- 
paigns revealed that all X-ray flares seem to have a simultaneous 
NIR counterpart (Eckart et al.||2004| |2006a 2008 a; Hornstein 
et al. 2007; Yusef-Zadeh et al. 2006a 2009, Dodds-Eden et al. 



Table 1. Observation log. 



2009, 1, while not all NIR flares have detectable X-ray counter- 



parts ( [Hornstein et al.||2007[ ). Considering also the hnear po 
larization properties of NIR flares (e.g., [Eckart et al.||2006b| , a 
general consensus emerged in which the origin of NIR flares 
was ascribed to synchrotron radiation of transiently accelerated 
electrons and X-ray flares to the corresponding inverse Compton 
emission, especially synchrotron self-Compton (SSC) (e.g., Liu 



[et a l. 2006a b). However, the recent study of an exceptionally 
bright X-ray/NI R flare detected on Apri l 4'^, 2007, chall enged 
the SSC model ( [Dodds-Eden et al.|2009|[Trap et al.|20T0| i. This 
proved the necessity to measure contemporaneous spectra (with 
good statistics) in infrared and in X-rays, to constrain the over- 
all spectral energy distribution (SED) of the flares. Yet, past at- 
tempts to obtain NIR spectral information on Sgr A* have given 
various and sometimes conflicting results (e.g., Gillessen et al. 
[20061 [Hornstein et al.|2007] l. 

On the other hand, broadband campaigns found indica- 
tions that contemporaneous X-ray/NIR flares could be accom- 



panied by delayed radio activity (Yusef-Zadeh et al. 2006a 
2006a| l. This was interpreted as a sign of expand- 



Eckart et al. 



ing plasmoids in the vicinity of the supermassive black hole, 
releasing optically thin and thick synchrotron emission whose 
peak evolves with frequen cy ([van der Laan 1966 1. Recently, 
Marrone et aL[ ( [2008) 1 and [Yusef-Zadeh et al.| ( [2008] l clearly 
detected an increase in the mm and sub-mm emission, using 
the Sub-Millimeter Array (SMA) and the Caltech Sub-millimeter 
Observatory (CSO) respectively, following an X-ray/NIR flare 
detected by Chandra and the Keck observatories. In the frame 
of an SSC/expanding blob model, Marrone et al. ( 2008| l argued 
that there should be an observable correlation between the X-ray 
to NIR peak ratio and the sub-mm time delay, which could help 
test this picture. 

We present here a new X-ray, NIR, and sub-mm cam- 
paign involving the X-ray Multi-Miror (XMM)-Newton satellite, 
the Very Large Telescope (VLT), and the Atacama Pathfinder 
Experiment (APEX) observatory, respectively. The data obtained 
enable us to investigate the aforementioned issues by measuring 
Sgr A*'s time variability and spectrum across a broad range of 
wavelengths. Concurrent public GeV y-ray data from the Fermi 
satellite are also examined. Additional radio data will be re- 
ported elsewhere. 

This paper is organized as follows. In Sect. |2] we give brief 
summaries of the data acquisition and reduction of each tele- 
scope. In Sect. |3] we describe the respective light curves of 
Sgr A*, identify new flares, and compare their properties to past 
detections. Sect. [4) is dedicated to the NIR color measurements 
performed on one particular flare, while Sect. |5] discusses pos- 
sible phenomenological interpretations of the light curves and 
SEDs of all the newly detected flares. 



2. Observations and data analysis 

We observed the GC in early April 2009 in the context of a 
large multi-wavelength campaign that was initiated by an XMM- 
Newton/VLT joint program and followed by an APEX monitoring 
scheduled accordingly. A detailed observation journal for this set 
of instruments is given in Tab. [T] Fermi data too were recorded 
at the epoch of this campaign as part of the ongoing monitoring 
of the GeV y-ray sky. We present these data and their respective 
treatments below. 



Orbit 


ObsID 


Start date 
[UTC] 


Exp. 
[ks] 


Band 


XMM-Newton/PN 


1705 
1706 
1707 


554750401 
554750501 
554750601 


2009 Apr 1,01:18:30 
2009 Apr 3, 01:56:31 
2009 Apr 5, 03:53:34 


37.8 
42.0 
32.4 


2-10 keV 
2-10 keV 
2-10 keV 


VLT/NACO 




082.B-0952A 
082.B-0952A 
082.B-0952A 
082.B-0952A 


2009 Apr 1,06:46:59 
2009 Apr 3, 07:25:41 
2009 Apr 3, 07:28:42 
2009 Apr 5, 08:44:10 


11.5 
9.7 
9.7 
6.6 


3.8 yum (L') 
2.1pm (K,) 
3.8 fim (L') 
3.8 fim (L') 


APEXfLABOCA 




183.B-0100A 
183.B-0100A 
183.B-0100A 


2009 Apr 1,05:31:43 
2009 Apr 3, 11:07:43 
2009 Apr 5, 06:40:10 


30.7 
8.2 
25.6 


870 pm 
870 pm 
870 pm 



2.1. XMM-Newton (2-10 keV) 



The X-ray satellite XMM-Newton Pansen et al.[[2001| l pointed 
toward the GC on three occasions between April and Apr 
5* 2009 (~3 X 40 ks). In the following, we will consider the 
data recorded with the Medium filter by the European Photon 
Imaging Camera (EPI C), i.e. flie PN ([Striider et al.||2001| l and 



the MOS 1-2 cameras ( [Turner et al. 2001 1 instruments, all op- 
erated in Prime Full Window mode. The data were reduced 
and analyzed with the Science Analysis Software (SAS) pack- 
age (version 8.0.0) and the latest calibration files available at the 
time of the analysis. They were not contaminated by soft proton 
flares from solar activity. Starting with the original Observation 
Data Files, event lists were generated through the emchain and 
epchain routines, and filtered with PATTERN <= 4 or 12 (for PN 
and MOS respectively), to ensure the selection of X-rays and 
rejection of cosmic rays. In order to prevent the use of X-ray 
photons that fell on bad pixels, we also applied the #XMMEA_EP 
and #XMMEA_EM filters on the event lists. Finally, artefacts from 
the calibrated, and concatenated datasets and events near CCD 
gaps were discarded by setting FLAG as selection criterion. 
We note that the high interstellar absorption toward the GC 
(A^H ~ 10^^ cm"^) prevents the use of photons below ~2 keV, 
so that the following count rates, fluxes, and luminosities will 
refer to the 2-10 keV band. 

The image obtained by the PN detector and reproduced 
in Fig. [T] (middle panel), reveals the usual constituents of the 
Galactic nucleus visible via XMM-Newton: extended thermal 
emission from the shell of the supernova remnant Sgr A East 
(Sakano et al. 2004 | l and point-like emission conincident with 
the position of Sgr A*. In addition, these observations allowed 
to locate another point source at only ~27" from Sgr A*, which 
we identify with the weak transient CXOGC J174538.0-290022 
(~2 X lO^'' ergs"' at 8 kpc, ,Ponti et al.|2009 1. 

To search for time variability of Sgr A*, we first fitted the 
positions of the point sources by running the edetect_chain 
procedure. A light curve of Sgr A* was then extracted from a 
10" radius circular region centered on the fitted position in each 
EPIC detectors. These light curves were background subtracted, 
with the background regions selected following the prescriptions 
given in Kirsch[ ( [2005| l and the counts scaled to the same ex- 
traction area as Sgr A . Finally, the 2-10 keV background sub- 
tracted PN, MOS 1, and MOS 2 light curves were summed and 
rebinned with a sampling of 600 s for the I'*' and 3"* observa- 
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Fig. 1. Multi-wavelength images of the Galactic center in early April 2009. North is to the top, East to the left, and the Galactic 
plane runs from upper left to bottom right. The sub-mm and X-ray images {left and middle panels) are summed over one or several 
observations, whereas the L'-band NIR image (right panel) corresponds to the summit of Flare A. Color scales are given in units of 
Jybeam counts pixel ' , and ADU, from left to right, respectively. 



tions, and 400 s for the 2"^. Note that, as resolved by previous 
Chandra deep observations, the persistent X-ray emission from 
a 10" radius region centered on SgrA*, not only contains the 
X-ray counterpart of SgrA*, but also includes persistent com- 
pact sources such as the star group IRS 13 (Baganoff et al. 2003) 
and the pulsar wind nebula candidate G 3 59.95-0.04 (,Wang et aL] 
|2006| l, possible transient point sources ( Muno et al.||2005| l, as 
well as a diffuse component ( Muno et al. 2003| l. 



2.2. VLT (2.1 fim and 3.8 ijm) 

The three X-ray observations described in the last paragraph 
were partly covered by NIR observations with the NACO instru- 
ment, NAOS (Nasmyth Adaptive Optics System, |Lenzen et al. 



2003 I + CONICA (COude Near Infrared CAmera, 



Rousset et al. 



2003 I, mounted on the fourth unit telescope, Yepun, at ESO's 
VLT in Paranal, Chile. These observations were divided into 
three half nights (April P', 3"^, and 5* 2009), which all started 
in the middle of the night and were extended as long as possible 
until sunrise. 

We acquired data solely in the L'-band (3.8 fim) during the 
and 3'^'' nights, while during the 2"*^ night, atmospheric con- 
ditions were favorable to also collect images in the /T-short (Ks) 
band at 2.1 jjm. During this latter night, we attempted a novel fil- 
ter switching technique with NACO in order to study the color of 
SgrA* from broadband quasi-simultaneous NIR measurements. 
Up to now, it has only been possible to investigate the spec- 
trum of SgrA* at the VLL through the integral field spectrograph 
SINFONI over relatively narrow bands: the H+K-hand (~1.7- 
2.4 jum, |Eisenhauer et al.||2005| l and the ^T-band (~2. 0-2.4 fim, 
[Gillessen et al.||2006| l. Using NACO, we obtained for the first 
time spectral data over a wider wavelength range and with better 
Strehl ratios, by constantly cycling back and forth between 
and L' photometric filters throughout the night. We carried out 
successive full switches (switching K^- to L'-band, acquiring two 
L'-band images, switching L'- to /Ts-band and acquiring two Ks- 
band images) on a timescale of only ~5-5.5 min. All the data ex- 
amined here were collected with the supergiant IRS 7 as natural 
guide star and the infrared wavefront sensor to close the adaptive 
optics loop of NAOS. Images were individually sky-subtracted. 



flat-fielded, and corrected for dead/hot pixels (see Fig. [T] right 
panel). 

In Fig. [3] [4j and [s] we display the light curves of SgrA* 
corrected for extinction with the usual values A^^ - 2.8 mag 
and Al> = 1.8 mag, in order to allow comparison with previ- 
ous works (e.g., [Genzel et aLl[2003| |Dodds-Eden et al.|[2009] l. 
They were obtained via aperture photometry and caUbrated with 
the fluxes of nearby stars of known and stable brightness: all 
fluxes were scaled so that the median flux of the relatively iso- 
lated star S 65 always corresponds to an apparent Ks magni- 
tude of 13.7 mag d Gillessen et al.|2009| l and the median flux of 
IRS16NW to an apparent L' magnitude of 8.38 mag ( [Viehmarm 
et al. 2005). More details on the 2"'* night light curves are given 
in Sect. [4] The duration of the 3"' night light curve is reduced 
compared to the other two due to a technical problem with 
the telescope. We stress that, at the time of these observations, 
SgrA* was confused with the star S 17, whose contribution to 
the light curves is expected to be a constant value of ~5.4 mJy 
and ~3.5 mJy in the K^- and L'-band, respectively (see S 17 mag- 
nitudes reported in Ghez et aL||2005 corrected with the afore- 
mentioned values of A^). 



2.3. APEX (870 ^m) 

In parallel to the X-ray and NIR surveys, a total of ~20 hr of 
sub-mm GC data were taken at 870 fim (345 GHz). We used the 
Large APEX bolometer cam era (LABOCA , [Siringo et al.|2009| l, 
on the APEX telescope (Giist en et al.|2006 )l at Llano Chajnantor, 
Chile. LABOCA is an array of 295 composite bolometers op- 
erating in the 345 GHz transmission window. The angular res- 
olution (FWHM) of each beam is ~19" for a bandwidth of 
~60 GHz. LABOCA's layout leads to a double beam spaced dis- 
tribution of the individual beams in a hexagonal configuration 
over a 1 1.4' field of view (FOV). 

The observations performed on the P' and the 5* of April 
spanned ~8.5 and ~7.5 hr, respectively, with an interruption of 
1 hr and 15 min between UT 8:45 and 10:00. During this time 
interval, the GC passed above the telescope's elevation limit of 
~80°. Observations on April 3"^ only covered the period after 
this GC transit (UT 1 1 :00 - 13:30). Note that the APEX and VLT 
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sites have basically the same longitude, but, contrary to the VLT, 
APEX could continue observing past sunrise, till the elevation of 
the GC was too low to allow for reliable monitoring purposes 
(elevation between 30° and 40°, depending on weather condi- 
tions). The precipitable water vapor during the three runs was 
~2, ~0.9, and 0.2-0.4 mm, respectively. 

The GC was observed using "on the fly" maps with scan- 
ning angles of -15°, 0°, or -1-15° relative to the orthogonal di- 
rection across the Galactic plane. The maps were built with a 
scanning speed of 2 arcmin s ' and steps orthogonal to the scan- 
ning direction of 30". The total integration time of each map was 
315 s. Each GC map was followed by observations of G 10.62 
and G5.89, standard LABOCA secondary calibrators. The at- 
mospheric zenith opacity was determined via skydips every 1- 
2 hr. With this observing setup, Sgr A* flux was monitored every 
~8 min. 

The data were reduced using the BoA software package. 
Reduction steps on the bolometer time series include temper- 
ature drift correction based on two "blind" bolometers, flat- 
fielding, calibration, opacity correction, correlated noise removal 
on the full array as well as on groups of bolometers related by 
the wiring and in the electronics, flagging of bad bolometers, and 
finally de-spiking. Each reduced scan was then gridded into a 
spatial intensity and weighting map. These reduction steps were 
applied to the GC and both calibrators scans. Flux calibration 
was achieved using G 10.62, which has a known flux density 
of 33.4 ± 2.0 Jy for LABOCA (Siringo e t al.|2009| . The abso- 
lute calibration accuracy is roughly 10%. To infer the relative 
calibration error, we applied the calibration curve of G 10.62 to 
G5.89 and determined the dispersion of the measured G5.89 
fluxes over the observing period. This yielded a relative calibra- 
tion accuracy of ~3%. 

For the determination of Sgr A*'s light curve we first gener- 
ated a model of the sub-mm emission of the region by co-adding 
all calibrated and pointing drift corrected GC maps. In this high 
signal to noise map (see Fig.[T] left panel), we note the presence 
of several molecular clouds in accordance with previous sub- 
mm GC surveys (I Pierce-Price et al.|2000| l, and the point source 
Sgr A*, which was fitted with a gaussian and later subtracted 
from the co-added map. We finally re-reduced all GC scans and 
subtracted the model signal from the time serie of each bolome- 
ter The resulting maps only contain the point source Sgr A*. The 
light curves were eventually constructed by fitting a gaussian to 
each individual map (see Fig.[3j|4] and|5]). The error bars were 
estimated through the RMS of the calibrator G5.89 (-2.7 %). 
Note that Sgr A* does not suffer from any significant interstellar 
extinction at 870 jim. 

2.4. Fermi (0.1-200 GeV) 

The Fermi space observatory carries the Large Area Telescope 
(LAT, [Atwood et al.|2009] l, a high-energy y-ray pair-conversion 
telescope covering the ~30 MeV to ~300 GeV band with a wide 
FOV (2.4 sr). It operates almost as an all-sky monitor and can 
provide regular ~50 min uninterrupted exposures on the GC. 
Using the public Fermi data archive, we explored the periods 
covered by the aforementioned observations of April 2009, and 
searched for a potential y-ray flaring activity correlated with 
variability in other spectral domains. 

The LAT data were analyzed with the Fermi Science Tools 
and following the analysis threads distributed by the Fermi 
Science Support Center We created events cube with the 
gtselect routine by selecting only y-ray photons ("diffuse" 
class events, N°3). We then used gtmktime and the "spacecraft" 



Fermi /LAT 1-200 GeV 




10 20 30 40 50 

Fig. 2. Map of the Galactic center built with 1-200 GeV y-ray 
photons collected by the LAT instrument during the P' year of 
the Fermi mission (Aug 2008/Aug 2009). North is to the top. 
East to the left, and the Galactic plane runs from upper left to 
bottom right. Each pixel has a size of 0.1° and the image was 
smoothed with a gaussian kernel of 0.15° radius. The color bar 
is expressed in counts and the positions of the point sources were 
taken from |Abdo etaL] ( |20T0l l. 



file which keeps track of the satellite's attitude as a function of 
time, to perform good time intervals selections. A zenith-angle 
cut was thereby applied (Zmax = 105°), rejecting y-rays emanat- 
ing from Earth albedo. The light curves were derived via aper- 
ture photometry (gtbin) with a circular extraction region of ra- 
dius 1°, centered on the position of Sgr A*. The fluxes thus mea- 
sured include the diffuse y-ray emission associated with cosmic 
ray interactions in GC molecular clouds. The number of counts 
in each bin, n, being governed by Poissonian statistics, Icr (84% 
confidence level) error bars were approximated by 1 -H V« + 0.75 
( Gehrels| 1986). Finally, we computed the exposures of each time 
bin in s~^ cm"^ with gtexposure. The resulting 0.1-200 GeV 
light curves are presented in the top panels of Fig. |3] |4] and [5] 
with a time bin of 600 s (10 min). 

For illustration purposes, we produced an average LAT im- 
age of the GC (Fig. |2f using all the data from the P' year of 
full sky survey, from August 2008 to August 2009. Only high 
energy photons from the 1-200 GeV band were considered to 
optimize the single photon angular resolution (< 0.6° at ener- 
gies > 1 GeV). 



3. Light curves and flares 

In this Section, we discuss the properties of the multi- 
wavelength light curves and compare them with the relevant Ut- 
erature. 
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Fig. 3. Light curves of SgrA* and its associated backgrounds 
(see text for details) during the P' observation. 



3.1. NIR 

The P' night of observations was marked by a strong NIR flare 
detected in the L'-band between 440-515 min past midnight 
(Fig. |3]l, hereafter referred to as "Flare A". It was brighter than 
the star S 2 (see the pale red dot north of SgrA* in Fig. [T] right 
panel) and peaked at 19 + 1 mJy (background subtracted), which 
makes it the 2"'^ brightest L' flare covered in X-rays and one of 
the few powerful NIR flares so far detected that are on the order 
of 20 mJy (in all NIR bands considere d) or beyond (|Dodds-Eden 
et al |2009[[Eckart et al.|200 8a b ; Mey er et al.|2006b||2007[ ). Itis 
worth noting that Flare A has a complex non-symmetrical mor- 
phology, characterized by a very sharp rise (a continuous gain of 
~19 mJy in ~7 min) followed by a slower and non-monotonic 
decay. This rise time is one of the fastest ever reported and is 
actually rather unusual among strong NIR flares. It hence places 
a stringent upper limit on the size of the overall flare emitting re- 
gion of ~lQRs (or ~1 AU), with a causality argument (R < c At) 
and considering a mass of 4 x 10^ Mq for the black hole. 

During the 2"'' night, the time resolution of the light curves 
was not as good as during the other nights, given the continu- 
ous cycling through the K and L filters, and hence makes the 
identification of flares more difficult. Sgr A* seems to have been 
continuously variable in both bands at a significant level (see 
Fig. |7] top panel, and § |4|, at least during the first two thirds 
of the observation. Our X-ray coverage allows us to clearly iso- 
late a new eruption — Flare B — between ~520 and ~560 min af- 



ter 2009-04-03 00:00:00 UTC though (see pai-agraph (3^ and 
Sect. |4|l. This time, the shape of the flare's light curve is more 
symmetrical, with no time delay between the peaks of the K^- 
and L' -bands larger than ~5 min, the time resolution of the light 
curve. The peak fluxes of Flare B are 3.7 + 0.5 and 8.2 + 0.5 mJy 
(background subtracted) in the Kf^- and L'-band, respectively. 
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Fig. 4. Same as Fig.^for the 2"'* observation. 



No significant flare was recorded during the 3"* night, which 
had an observed mean flux density of ~4 mJy. 

3.2. X-rays 

No prominent X-ray flare (> 10 times the quiescent level) was 
detected in April 2009. We concentrate, hereafter, on the signif- 
icant peaks (>3cr) of the X-ray light curves detected in parallel 
to the NIR flares mentioned above. 

Flare A coincided with a very little increase of the XMM- 
NewtonfEPlC cameras count rate (see Fig. |3] top and middle 
panels). Though barely significant (3.5cr), this could be the X- 
ray counterpart to Flare A and thus would be the faintest X-ray 
flare so far identified with XMM-Newton, all previous ones be- 
ing, at least, a factor 2 brigh ter ([Goldwurm et al . 2003 ; Porquet 
et al. 2003; Belan ger et"aL]|2005[ [Porquet et al.||2008| l. When 
considering the overall flare, the X-ray excess occurs in the mid- 
dle of the NIR flare, which is consistent with the flares that have 



been entirely covered in NIR and X-rays over the past (|Eckart 
et al. 2006a 2008al [Dodds-Eden et al.|2009l l. In view of the low 
statistics of this X-ray signal, we cannot extract proper spectral 
information. Assuming a typical photon spectral index, F = 2 
(with dN/dE oc E'^ phs"' cm^^keV"'), and a typical column 
density, A^h = lO^-' cm"^, we estimate through WebPIMMS that 
Flare A reached an unabsorbed, background subtractecQ flux 
of 2.2 X 10"'^ ergs"' cm"^ in the 2-10 keV band. This trans- 
lates to a peak luminosity of 1.7 x lO"*"* erg s~ ' at 8 kpc, i.e 



-7 times the quiscent level (2.3 X 10^^ ergs ' , |Baganoff et al 



2003 ). This very low significance X-ray flare associated with a 
relatively bright NIR flare is quite surprising in regard to the 



' An average level of quiescent count rate including the constant con- 
tibution of neighboring sources (§ |2.1[ l, was removed as in |Goldwurm| 
etaDpOOSll. 



6 G. Trap et al.: Concurrent X-ray, NIR, Sub-mm, and GeV y-ray observations of Sgr A* 

Table 2. Properties of the X-ray/NIR flares observed so far. 



References!^ 


Date 


Instruments 


NIK FILTER 


ix/i[|uics|^ 


fNIR[] 
[mJy] 


^nirU 
[mJy] 


Id. 


[1] 


1 Apr 2009 


XMM 1 VLT 


3.8 /jm (L') 


< 7 


19 




A 


[1] 


3 Apr 2009 


XMM 1 VLT 


3.8 //m (L') 


10 


8 




B 


[2,3] 


4 Apr 2007 


XMM 1 VLT 


3.8 (L') 


103 


26 




C 


[2,4] 


4 Apr 2007 


XMM 1 HST 


1.7 fim 


26 


7 


16 


D 


[2,4] 


4 Apr 2007 


XMM/HST 


1.7 //m 


37 


4 


9 


E 


[2,4] 


2 Apr 2007 


XMM /HST 


1.7 /jm 


14 


10 


22 


F 


[5,6] 


31 Aug 2004 


XMM/HST 


1.6 jum 


< 40 


10 


24 


G 


[7] 


31 Jul 2005 


Chandra / Keck 


3.8 //m (L') 


< 1.2 


15 




H 


[7] 


17 Jul 2006 


Chandra / Keck 


3.8 Aim (V) 


20 


> 5 




I 


[8] 


30 Jun 2005 


Chandra / VLT 


2.1 /jm (ATs) 


3 


14 


4 


J 


[9] 


7 Jul 2005 


Chandra / VLT 


2.1 //m (K,) 


15 


11 


3 


K 


[10] 


19 Jun 2003 


Chandra / VLT 


2.1 /jm (/fs) 


3 


> 7 


> 2 


L 



[6] 

b 



References: [1] This worlc, [2] 
i 2006a| , [7] 



Yusef-Zadeh et al. 
The //5r 1.6 and 



Porquet et al. ( 2008b, [3] 



Dodds-Eden et al. 



Hornstem et a . 

jands correspond to the NICMOS filters F160W and F 170M, respectively. 



|2009| , [ 4] I Yusef-Zadeh et al | ( |2009|, [5]|Belanger et al.| ( |2005 



1 [2007| , [8]'Eckar t et al.| l |2d68a^ , [9] |Eckart et al.| ( |2006a^ , [10] 



Eckart et al. 



ergs 



Baganofretal.|2003| 



7 jim pass 

Peak X-ray luminosity (2-10 keV) in units of the quiescent level (~2 x 10^' 
when the peak flux was missed or not simultaneous to the NIR peak. 

Background subtracted peak NIR fluxes in the bands indicated in the fourth column, dereddened assuming Ai 
5.03 mag, A2.i;jn, = 2.8 mag, and A3 g^,^ = 1.8 mag. 
' Equivalent NIR flux in the L'-band extrapolated assuming a power-law spectral index a = - 1 (see Sect.|4]l. 



1 2004 1, 



>' and '<' signs were added 
4.5 mag, A, = 



° 

i -2 
~ 0.3 

g 0.25 
o 

6 0-15 
0.1 



a- 



- Fermi; LAI, 0.1-200 GeV 



- XMM / EPIC, 2- 




- VLT / NACO, 3.8lim 



APEX / LABOCA, 870 urn 



300 400 500 600 700 

Time (min from 2009-04-05 00:00:00 UTO) 



Fig. 5. Same as Fig. [3]for the 3'^'^ observation. 



bright X-ray flare (~100 x the quiescence, Porquet et al.|2008] l 
detected simultaneously to a strong L' NIR flare (~30 mJy) 
in April 2007 (Dodds-Ede n et al.||2009[ ). It is however some- 
what reminescent of the X-ray non-detection by the Chandra 
observatory of a ~15 mJy L' flare observed in July 2005 by 
|Hornstein et al.| ( [2007[ ). Yet, we stress that the point spread func- 
tions (PSFs) and spectral responses of the Chandra/ ACIS and 
XMM-NewtonfEPlC instruments differ, which prevents straight- 
forward comparisons. Actually, the fact that the EPIC's PSF is 



broader than ACIS's makes the XMM-Newton satellite less sen- 
sitive than Chandra to soft (F x 2) and faint X-ray flares. The 
background subtracted peak XMM-Newton count rate of Flare A 
converts, indeed, to a Chandra/ AClS-l rate of ~0.034 counts s"' 
in the 2-8 keV band. This would have lead to a clear Chandra 
detection, even more significant than two other faint X-ray/NIR 
flares observed with Chandra by |Eckart et al. (2004 2008a I. 

Flare B, on the other hand, was more significantly detected 
by XMM-Newton and observed to last longer (~40 min), in sync 
with the NIR light curve (see Fig. |4] top and middle panels). 
Using the same method as above, we estimate its background 
subtracted and unabsorbed peak flux is 3.2 x 10"'^ erg s"' cm"^ 
in the 2-10 keV band, i.e. 2.4 x 10^** ergs"' (-10 x the quies- 
cence). Note that, with Chandra, it would have been detected at 
a rate of 0.049 counts s"' (2-8 keV), which is still faint com- 
pared to the other X-ray flares with a NIR counterpart known as 
we discuss below. The fact that the NIR flux was relatively high 
about 100 min before Flare B with no simultaneous detectable 
X-ray emission is striking and could be another example of a 
"NIR only flare" pornstein et al.||2007| l or indicate that NIR 
flares can start earlier than the X-ray ones ( |Dodds-Eden et aL] 
i2009| ). 

Comparing Flare A and B, it is manifest that the ratio of 
the peak NIR flux over the peak X-ray flux is highly variable 
from one flare to another. It is known that the NIR flares oc- 
curence rate is higher than the X-ray one, so it was expected 
that there should exist NIR flares with very little or even no de- 
tectable X-ray emission. But a somewhat counterintuitive fact 
that emerges now is that the brightest NIR flares are not nec- 
essarily accompanied by bright X-ray flares. Indications of that 
result had been obtained with the analysis of the July 2005 and 
July 2006 flares in Hornstein et al. (2007 1, although the peak 
NIR flux of the July 2005 event was missed. Similarly, Yusef-| 
Zadeh et al. ( 2009[ ) showed in their Fig. 21 (a) and (b), with a 
consistent XMM-Newton/HST dataset, that the NIR to X-ray ra- 
tio fluctuates. 

To allow better comparisons between flares concurrently ob- 
served in X-rays and NIR, and visualize their disparities, we 
give an inventory of all the results published thus far in Tab. |2] 
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10' 



10' - 



10" 



X-ray quiescent state 



10' 



Fig. 6. Relation between the X-ray and NIR peak fluxes of afl 
the flares so far covered with muhi- wavelength observations. See 
Tab.|2]for details about each individual flare. Events A and B are 
the ones reported in this work. X-ray upper limits and NIR lower 
limits are indicated by down pointing and right pointing arrows, 
respectively. 



and plot the X-ray peak fluxes Vs. peak NIR fluxes in Fig. [6] 
Despite a relatively easy comparison between the measurements 
conducted by XMM-Newton and Chandra, NIR comparisons be- 
tween various authors are more diflicult because of different 
observational setups and flux corrections (e.g., background re- 
moval, extinction correction). Flares A, B, C (Tab. [2| are still 
easy to compare, since they were observed by the same in- 
struments in the same pass-bands, at roughly the same epoch. 
Flares H and I can also be added to this group, provided a 
background correction is applied and consistent extinctions em- 
ployecj^ Finally, peak NIR flux of the remaining seven flares (D, 
E, F, G, J, K, L) need to be extrapolated to the L'-band, assum- 
ing a particular infrared spectral index (see Sect. |4]l, to ensure 
consistency with the other flares. 

As shown in Fig. |6] no particular correlation between the X- 
ray and NIR strengths is apparent with the currently attainable 
statistics. Though the logarithmic scales of the plot axes do not 
reveal it. Flare C really stands out of the group as a strong event 
in both X-rays and NIR. 



3.3. Sub-mm 



On April P' 2009, sub-mm observations overlapped with the 
NIR and X-ray coverages (Fig. [3]l. Sgr A* is found to be con- 
tinuously variable at 870 //m, though on a longer time scale than 
at shorter wavelength^ We first note that there are no abrupt 



^ We used a L' background of 3.6 mjy as given in caption of Fig. 5 
and converted the extinctions quoted in caption of Tab. 1 ( |Homsteiri| 



et al. 2007 1. 

" Analogous variability, on timescales ranging from hours to days, 
was also reported at wavelengths longward of ~1 mm (see e.g., Wright 
& Backer|1993[|Bower et al |2002||Zhao et al.|2003l|Hermstein et al 



2004 



2007 



[Miyazaki et al.|2004[ [Mauerhan et al.|26o5| |Yusef-Zadeh et"aL 



changes in the APEX light curve at the time the VLT detects the 
rise of Flare A. This means that a simultaneous sub-mm coun- 
terpart to this NIR flare is not detected as such, and corroborates 
similar pas t non-detections ([Marrone et al.||2008[ |Yusef-Zadeh| 
et al.|2008t|Eckart et al.|20'08bl [ 

679+4 min after 2009-04-01 00:00:00 (UTC), the sub-mm 
flux of Sgr A* peaked at 5.4+0.2 Jy, and marks the summit of a 
dramatic sub-mm flare, among the brightest ever detected, that 
lasted for ~150 min and that we will refer to as Flare A'. When 
we subtract the background emission (taken as the minimum flux 
of the observation, 2.2 Jy), the peak is reduced to 3.2+0.2 mJy. 
The rise time of Flare A' is on the order of 20 min, which lim- 
its the size of the sub-mm flaring region to only ~3QRs with 
the same argument as the one used in § 3.1 Yet, this size con- 



straint will probably be revised soon when Very Long Baseline 
Interferometry (VLB I) measurements is achieved in the sub-mm 
range. Recent size measurements at 1.3 mm have indeed already 
yielded an intriguing upper limit to the quiescent size of Sgr A* 
of -40 juas ~4Rs (Doe leman et al.|2008| l. 

Interestingly, Flare A' lags Flare A by ~200 min and is not 
associated with any other flaring episode, at least as far as we 
can tell from the X-ray light curve. As summarized in Tab. [3] this 
new flare detected with high confidence adds to a handful of past 
detections (sometimes yet debatable) of short mm/sub-mm ac- 
tivity in the aftermath of X-ray and/or NIR flares. One way often 
proposed to relate these phenomena is that X-ray /NIR flares are 
associated with a plasmoid expansion radiating via synchrotron 



emission (van der Laan||I966 see Paragraph 5.1 1. If this plas 



moid is also responsible for the production of X-rays through 
SSC emission, then as pointed out by Marrone et al. ( 2008 ), there 
is a tight relation connecting the time lag between the X-ray/NIR 
peak fluxes (Fx, ^^nir) and the peak flux at a smaller frequency 
V. This lag is given by: 



„(p+2)/2 (4+p)/2 



i/a(3+2p) 



(1) 



where Bq is the inital magnetic field, p the spectral index of the 
radiating particles [dN/dE oc E"''), and 5 is a parameter de- 
picting the speed at which the radius R of the plasmoid grows 
(R oc f*). So far, only three flares were fully covered by contem- 
poraneous X-ray, NIR, and sub-mm observations and allow for 
a test of this relation: Flares A, H, and I. The comparison be- 
tween H and I was done by Marrone et al.'('2008), who showed 
a consistency between the data and the plasmoid/SSC picture. 
Yet, when considering Flare A/ A', this tendency of "the higher 
the ratio /^x/^nir> the longer the delay" no longer holds. Indeed, 
Flare I has a ratio about 10 times the ratio of Flare A, and still 
presents a delay about twice shorter This tends to indicate that 
X-ray/NIR flare and sub-mm flare emissions may not originate 
from the same particles (see also § |5.1| l. 

Regarding the other two APEX observing runs of April 2009, 
no sub-mm flares were detected, notwithstanding an obvious 
flux decrease observed during the last run. The minimum flux 
reached was 2.2 Jy, justifying our previous choice for the steady 
quiescent background level. 



3.4. GeVj-rays 

The Fermi light curves in the top panels of Fig.[3]|4] and|5] are 
characterized by frequent data gaps on the order of 2 hr during 
which the GC is outside of the LAT's FOV The statistics within 
each bin is remarkably poor; many 10 min bins do not contain 
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Table 3. Short mm and sub-mm flare delays relative to X-ray or 
NIR flares. 



References|^ 


Date 


Instruments 


Time lag|^ 
[min] 


[1] 


1 Apr 2009 


APEX (870 /jm) 


200 


[2,3,4] 


17 Jul 2006 


SMA (1.3 mm) / CSO (850 /jm) 


100 


[2] 


31 Jul 2005 


SM4(1.3mm) 


20 


[5] 


3 Jun 2008 


APEX (870 ^^m) 


90 


[6,7] 


7 Jul 2004 


SMA (890 /im) 


< 120 


[4] 


5 Apr 2007 


SMA (1.3 mm) 


160 


[8] 


3 Sep 2004 


CSO (850 /jm) 


160 



References: [1] This work, [2] 



Marrone et al. 



[Zadeh et al.' ('2008'), [4] Yuse f-Zadeh et al.| j2009b, [5] pclcart eT^ 



(,2008bl, [ 6] Eciiart et al. (2006a|l, [7] |Eckart et al 
Zadeh etal.J|2006ai. 



12008 



p009 



[3] 



Yusef- 



[8] lYusef- 



Time by which the sub-mm peak lags the X-ray or NIR peak. 



£ 4 



VLT/ NACO, 3.8nm (L) 
VLT / NACO, 2. 1 2(im (K) 




480 500 520 540 560 580 
Time (min from 2009-04-03 00:00:00 UTC) 



any photon at all. No significant peak in the light curves is de- 
tected. 

Though Flare A and most of Flare A' fell outside the cover- 
age of the GC by the LAT, the rising flank of Flare B was coin- 
cident with the end of a LAT exposure. Zero photons were de- 
tected in the time bin at the peak of Flare B. Given the Poissonian 
statistics, the 3cr (99.87% confidence level) upper limit on the 
count rate of any source in the extraction region is 6.6 photons 
( |Gehrels|1986| ) for the bin considered, which has an exposure of 
6.2 X lO'* s~^ cm~^. This converts to a flux upper limit for Flare B 
of 1.1x10-4 phs ' 



nosity upper limit of 5.5 x 10"" erg s 



'cm--, or 7.2x10" 

38 



erg s cm , and a lumi- 
(0.1-200 GeV)at8kpc. 
This is not very constraining in regards of the flux 
quoted for the GC Fermi centralpoint source, 1FGLJ1745.6- 
2900, in the P' year catalogue^] ( Abdo et al.||2010| : 1.7 x 
10 phs 'cm^^ (0.1-100 GeV) for a photon spectral in- 
dex of ~2.2. Fermi thus confirms the presence of a GeV 
source in the GC, as once detected by the Energetic Gamma- 
Ray ExperimenT (EGRET) onboard the Compton Gamma-ray 
Observatory ( Mayer- Hasselwander et al. 19981. The origin of 
IFGL J 1745. 6-2900 is unknown. It is quite possible that this 
GeV emission is produced in GC molecular clouds just as the 
TeV emission of the same region (Ballantyne et al.|2007i Melja] 
& Fatuzzo 2010 i. But Sgr A* itself is also of course a valid 



candidate, given the positional coincidence of the LAT source 
and the supermassive black hole: IFGL J 1745. 6-2900 lies at 
the galactic coordinates [l = 359.941°, b = -0.051°] with a 
95% confidence radius of 1.1', which is only ~21" away from 
Sgr A*. One way to unambiguously attribute this source to the 
supermassive black hole could have been the detection of a flar- 
ing activity correlated with the X-ray one. However, we did not 
detect any GeV signal with the sensitivity we achieved dur- 
ing Flare B (which limited our detection capability to large 
GeV flares of amplitude at least ~65 times the steady flux of 
IFGL J 1745. 6-2900). The present study is hence analogous to 
the unfruitful previous attempts to discover correlated emission 
between the X-ray flares of Sgr A* and the flux of the y-ray 



sources I GR J 17456-2901 an d HESS J1745-290 (see [Trap et alT 
[2010 ; A haronianetal.|2008| for the INTEGRAL and HESS non- 
detections, respectively). 



These measurements were obtained with a modeling of the under- 
lying diffuse emission. A variability index of 18 was measured taking 
11 time segments into account and implies that IFGL J1745. 6-2900 is 
non-variable on a monthly timescale. 




2 3 4 

K-band fiux density (mjy) 



Fig. 7. {Top panel) NIR obervations of the 2"'' night background 
subtracted and corrected for extinction with Aii - 1.8 mag 
and Ak^ - Ai> + E{Kfi - L') - 3.3 mag. (Bottom panel) 
Corresponding spectral index a plotted against dereddened K^,- 
flux. Open squares indicate measurements made during Flare B. 



4. NIR colors of Sgr A* 

Here we describe the results of the spectral imaging of Sgr A* 
we made during the 2"'* night of VLTTNACO observations. 

To compute the NIR photon spectral index a, defined by 
Fy oc v", we produced the L'- and Ks-hwA light curves of Sgr A* 
presented in Fig. [7] (top panel), for which we averaged the pairs 
of images taken at each filter switch. The error bars incorporate 
the standard deviation of these values as well as the error de- 
termined from the standard deviation of a comparison star of 
similar flux. Data points were interpolated to the same times 
and a large aperture was used in both bands to fully include the 
contribution of S 17. A global background correction (including 
S 17) was made by subtracting the lowest flux value. It is im- 
portant to note that a is highly dependent on the excess color 
E(K^ - L') - Ak, — Al', since: 



log {\Q<^^':.-'^L')l2.i X F'j^^f/Flf) 



(2) 



G. Trap et al.: Concurrent X-ray, NIR, Sub-mm, and GeV y-ray observations of Sgr A* 
Table 4. Flare NIR spectral measurements. 



9 



Kefer e nce s I 1 


Dates 


— f 

Instruments 


g 

oANDS 




[1] 


3 Apr 2009 


VLT /NACO 


2.1, 3.8 fim 


-0.4 ±0.3 


[2] 


16, 31 Jul 2005, 2 May, 17 Jul 2006 


Keck /NIRCl 


1.6, 2.1,3.8,4.7 fim 


-0.6 ± 0.2 


[3] 


2, 5 Apr 2007 


HST /NICMOS 


1.45, 1.7 /im 


-0.7 ±0.2,-1 ±0.3+ 


[4,5] 


4 Apr 2007 


VLT/NACO, VISIR 


3.8, 11.8//m 


> -1 


[6] 


29 Apr 2005 


&cyt/ OSIRIS 


2.02-2.38 fim 


-2.6 ±0.9 


[7] 


18 Jun 2005 


VLT/ SINFONI 


2-2.45 /im 


-(0.6-1.3) ±0.2"^ 


[8] 


15, 17 Jul 2004 


VLT/ SINFONI 


1.7-2.45 //m 


-4± 1 


[9] 


26 Jul 2004 


Keck /mRC2 


2.1, 3.8 fim 


-0.5 ± 0.3 



References: [1] This work, [2] 



Hornstein et al. 



2007 I, [3] 



. ||Yusef-Zadehet al.|j2009|, [4]|Dodds-EdenerZl|2009t , [5] |Trapetai1|2010t , 

[6] |Krabbe et"aL] ( (2006| . [7]|Gillessen et al.| j2006f ," i8] |Eisenhauer et al.| ( |2005| ir[9] |Ghez et al.| l |2005[ >. 
' a is defined by Fy oc v".^ indicates when there is a correlation observed between the flux and the spectral index. In the case of Yusef-Zadeh 



et al. 1 2009b, only the brightest flares are reported. 



where F are the reddened fluxes. Using known early- type 
star^in the vicinity of Sgr A*, we determined an average color 
excess E{Ks-L') - 1.5±0.1 (see e.g. Hornstein et al. 2007| l. The 
corresponding values of a are plotted in Fig.[7](bottom panel) as 
a function of the /Ts-band flux. Given the statistics of the data, no 
particular correlation between a and the flux is found. The mean 
value obtained over the duration of Flare B is a = -0.4 + 0.3. 

For comparison, we list in Tab. [4] all the previous at- 
tempts to measure the spectral indices of flares from Sgr A*. 
The different values of a thus obtained are however diflicult 
to compare, since each of these works either probed different 
spectral ranges/flux amplitudes/epochs or used different instru- 
ments/techniques/reddening laws/backgrounds. Owing to the 
experimental setup used in this study, our result can only be eas- 
ily compared to the work of Hornstein et al. ( 2007| l. Likewise, 
these authors made use of a ground telescope with adaptive op- 
tics and the filter cycling technique, though they also accumu- 
lated data in the H- and Mj-bands, and used Keek's laser guide 
star. They found a constant "universal" value for the index of 
several flares, no matter the instantaneous flux, of a = -0.6±0.2, 
which is fully consistent with our result. Such a stable index 
seems to support a sort of a IR "quiescence" on average even 
though this is debated ( jPodds-Eden et al.|20I0a| l. 

It seems that all the values for a recently obtained suggest 
a > -1, which translates toyS > (with vFy oc v^), thus justify- 
ing that "infrared flares are blue" and arise from a particle pop- 
ulation distinct from the one producing the quiescent "sub-mm 
bump" ( Dodds-Eden et al.|2009 Trap et al.|2010 1. Presuming a 
synchrotron origin for the NIR flares and the radiating electrons 
distributed in a power-law of index p - \ - 2a (dN/dE oc E^''), 
this suggests that p < 3, which has important implications for 
the fit of flare spectra and in particular the SSC components (see 
Paragraph |5.2| l. 



5. Discussion 

Below, we discuss a simple modeling of Flare A/ A' light curves 
and SED based on the expansion of a plasmoid releasing syn- 
chrotron and SSC radiation. Flare B SED is then analyzed in 
the context of different non-thermal emission processes. Models 
for SEDs are essentially static and only meant to show in which 
direction our data may constrain the theory of Sgr A*'s flares. 



5 We used S 2, S 65, S 66, S 72, and S 87 which have K magnitudes 
between 13.5-14.6 mag and are spectroscopically identified to have 
spectral types between 09-B2 (Eisenhauer et al. 2005[ l. 



Evolutionary effects and the physical conditions leading to such 
scenarios will have to be studied in more details in future works. 



5.1. Plasmoid expansion 

Over the past years, several authors have employed the adiabatic 
expansion of a plasmoid (|van der Laan 1966) in the context of 
cm to sub-mm variability of Sgr A*, often observed to follow 
X-ray/ NIR flares (e.g., Yusef-Zadeh et al ||2006b| |Eckart et al.| 
|2006a| and the references listed in Tab. [3|). The use of this ex- 
pansion model for the flares of Sgr A* is motivated by: 

- the timelags separating X-ray/NIR flares from sub-mm flares 
(Tab. [3]); 

- the timelags separating X-ray/NIR flares from radio flares 
( Yusef-Zadeh et al.|2009t Kunneriafli et al.|2008) ; 

- the ti melags between different radio bands within the same 
flare ( |Y\isef-Zadeh et al.|2006bF 



- the increase of linear polarization as sub-mm flares reach 
their maximum ( |Marrone et a l. 20(Wy, 

- the similarity between the flare timescales in X-ray/NIR and 
mm/sub-mm, which suggests that energy losses are not dom- 
inated by radiative cooling, with expansion a good alterna- 
tive candidate^ 

Fig. [8] captures the essence of the model and gives the re- 
spective initial positions in the spectral domain of the wavebands 
used in this work. A ball of plasma of initial radius Rq, is filled 
with relavistic electrons of uniform volumic number density «e, 
arranged in a power-law (d«e/dr = k^y^P, ymin < y < Tmax)- 
This plasmon is pervaded by a magnetic field, B, and yields an 
initial synchrotron spectrum made of two connected power-laws 
(Fig. [8] middle panel, red curve). Below the initial synchrotron 
self-absorpti on (SSA) frequency vq (oc {Rok^}^'^i'+'^'B^''^^>^^P^*\ 
|Gould|1979| l at which the flux attains its maximal value, Fq, the 
spectrum is optically thick (Fy oc y^/^). For v > vq, the spectrum 
is optically thin (Fy oc v(i-/')/2) and breaks exponentially at a fre- 
quency proportional to y^^j^B. The dependence of the flux, Fy, 
on the radius, R, is given by: 



Fy = Fo 



5/2 



R\ 1 - exp(-Tv) 
RqJ l-exp(-To) 



(3) 



* Detailed models of these observations seem to exclude the sim- 
plest expansion scenario though, because the observed sizes of the radio 
emitting regions are bigger than van der Laan predicted sizes ( jFalcke| 
et al. 2009: Maitra et air2009| ). 

' Radiative cooling effects may however also affect the expansion 
model in the mm range i Li et al. 2009b. 
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(1966 1 model. As the plasmoid expands {left panel), its initial self-absorbed synchrotron 
spectrum is shifted towards smaller frequencies (middle panel). If the source is observed at optically thin frequencies from the 
beginning, then it remains optically thin and the flux keeps decreasing as time moves on (right panel, NIR curve). Alternatively, 
if the source is initially optically thick, then the flux rises until the source turns transparent and subsequently declines in the thin 
regime (right panel, sub-mm and radio curves). 



where the optical depth, Ty, obeys the equation: 

-{p+4)/2 , p v-(2/j+3) 



Tv = To 



'{- 



(4) 



and where the optical depth tq at vo is set by the condition that 
Fy be maximal [dy Fy(v{), t) - 0], or: 



exp(To)-To(pH-4)/5-l = 0. 



(5) 



In the following, we adopt a constant velocity, Vexp, for the ex- 
pansion so that: 



R(t) = Ro + VexpCf - to) . 



(6) 



To test this model with our multi-wavelength data of April 1, 
we propose two phenomenological scenarios (I and II) with dif- 
ferent sets of values for the parameters (p, Rq, vq, Fq, V^xp, k), 
see Tab. [5]). To assess the viability of these models not only with 
the light curves but also spectrally, we calculate the plasmoid 
initial synchrotron spectra with particular values of ymin, Tmax, 
B, and n^., and also consider the naturally associa ted SSC emis- 
sion, using the formalism of |Krawczynski et al.| (|2004p] In all 
cases, the NIR emission is assumed to be optically thin, and so 
^submm < vq < ^^NiR (Fig |8] middle panel). The NIR cou nterpart 
of Sgr A* is indeed known to be highly linearly polarized ( Eckart 



let al.|2006b[ [Meyer et al.|2006albl|2007[|Trippe et al.|2007 ') and, 
if it was optically thick, then an expansion model would anyhow 
predict a sub-mm peak flux lower than the NIR one (Fig. [8] right 
panel), which is inconsistent with the observations. 

Model I. The scenario that requires the least number of in- 
gredients to accomodate the NIR, sub-mm, and X-ray light data 
all at once, is a plasmon initially producing the NIR and X- 
ray flare by synchrotron and SSC, respectively. As it then ex- 
pands, the NIR/X-ray flare decays and the sub-mm flare arises. 
Given the great difference between the NIR and sub-mm fluxes, 
this model demands an extremely steep power-law distribution 
(p - 5) and a high maximum flux (Fq - 22 Jy) at the SSA 



In the progress of this work, we found a problem in the treatment 
of the SSA in H. Rrawczynski's code, which derived from a typo in 
the numerical expression in SI units of the absorption coefficient given 
bylLongair ( 19941, that should be written: Xv = 20.9 k' s(''+2)/2 (5.67 x 
10^)" b(p) v-^''^. 



turnover. We choose the different parameters so as to respect the 
flare size constraint (j?n < 10/?s) and the essential idea of the 
van der Laan ( 1966) model (v|^™ < vo, where v^™ oc y^^^Vg, 
= etl/\2nml.j is the gyration frequency, and nif. is the mass of 
the electron). Fig.|9](top and middle panels) shows that the over- 
all NIR/sub-mm timescales and flux characteristics of Flare A/ A' 
can be roughly reproduced (predictions at 100 //m and 3 mm are 
also plotted), and Fig. |9] (bottom panel) that the observed ini- 
tial SED constraints can be met. However, the initial steep NIR 
spectrum contradicts the measurements of a made with Flare B 
for example and exceeds the mid-infrared (MIR) upper limit ob- 
tained with VLrA'ISIR during the April 2007 campaign ( |TYap| 
et al. 2010) by a factor ~4. We therefore disfavor model I, al- 
though we cannot categorically exclude it given our lack of in- 
frared spectral information on Flare A. 

Model II. We now assume that Flare A and A' are not di- 
rectly connected by the same synchrotron power-law, and that 
Flare A simply acts as a trigger for a plasmoid expansion that 
will in term give rise to Flare A'. The plasmon power-law index 
can thus be chosen less steep (p - 2) and the couple of parame- 
ters (ymux, B) so that the initial plasmoid spectrum breaks in the 
far-infrared (FIR). It thus gives no significant synchrotron emis- 
sion in the NIR and X-ray bands, and the same holds for the SSC 
radiation. Yet, if the expansion starts at the beginning of the NIR 
flare, the relatively long time interval until the maximum of the 
sub-mm flare makes it difficult to obtain a narrow peak in the 
light curve at 870 yum. Therefore, we allow the time at which 
the expansion is initiated to be postponed by 140 min, which 
gives an excellent fit to the APEX light curve (Fig. [TOji. This is 
our favored model and it could be tested with radio millimeter 
data, since it predicts a distinctive 3 mm flare peaking at 0.6 Jy 
around 350 min after the onset of the NIR/X-ray flare. We stress 
that, though ad hoc at this stage, this delayed expansion can be 
physically motivated in a "coronal mass ejection" picture by the 
fact that it would take the accelerated particles ~200 min to be 
advected out of the inner parts of the accretion flow, where the 
high viscosity would prevent an immediate adiabatic expansion 
(see |Dodds-Eden et al.|20T0bl l. 

In all the above cases, we always used = ^s- But, as 
formulae |3] and |4] only involve the ratio R/Rq, we would have 
obtained the exact same sub-mm light curves by multiplying 
Rq by any factor x, provided the velocity Vexp quoted in Tab. |5] 
was multiplied by the same x. Since 10 Rs is an upper limit on 
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Fig. 9. (Top panel) APEX and VLT light curves plotted on the 
same scale. Theoretical predictions of model I for several wave- 
lengths are overplotted in blue. (Middle panel) Magnified view 
of the VLrdata. (Bottom panel) The synchrotron and SSC curves 
give the initial SED of the plasmoid in model I. They fit the sim- 
ulaneous NIR and X-ray data (in red and blue respectively) of 
Flare A. The non-simultaneous sub-mm data of Flare A' are also 
indicated in green for information. The references of the quies- 
cence data in light gray are given in Trap et al. ( 2010 1. 



the size of the flaring region, we find that the maximum expan- 
sion speed required is 0.05 c, which is consistent with the other 
sub-relativistic expansions found previously ( |Yusef-Zadeh et alT] 
2008||Eckart^eral . 2008b 2009). This means, incidentally, that 
the plasmoid cannot escape the gravitational potential well of the 
black hole, unless it has a relativistic bulk motion, inside of a jet 
for example d Falcke et al.|2009[partra et al.|2009) . 

The modeling we present here compares well with previous 



works involving similar spectral windows (Eckart et al. 



|2009| ), although we used only one plasmon compared to the 
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Fig. 10. Same as top panel of Fig. |9]for model n. 
Table 5. Parameters of Flare A and Flare B models. 



Flare A 



Flare B 



Parameters 


Model I 


Model II 


SSCl 


SSC2 


SB 


P 


5 


2 


2 


2.4 


2.4-3.4 


«o(«s) 


1 


1 


1 


1 


1 


vo(Hz) 


2x 10'2 


2x 10'- 








Fo (Jy) 


22 


17 








Vexp (C) 


0.0017 


0.005 








to (min) 


445 


588 








Tmin 


20 


1 


1 


1 


1 


Tmax 


2x 10^ 


100 


1 X 10'' 


2x 10" 


6x 10' 


B (G) 


90 


100 


2 


5 


21 


Jle (cm"'') 


1 X 10' 


3x 10' 


1 X 10' 


6x 10' 


5 X 10** 



multi-component models needed in these works. Nevertheless, 
our main result is that the simplest form of an abiabatic ex- 
pansion scenario (model I) is inadequate to reproduce all our 
NIR/sub-mm/X-ray data at once, so that we favor expansion 
models in which NlR/X-ray flares are not initiated by the par- 
ticles of the plasmon creating the sub-mm flare. This is in con- 
trast with the study of ,Eckart et al. ( |2008b_ ), in which models 
analogous to model I are generally favored. We note in particu- 
lar that "model BI" of Eckart et al.| ( |2008b"l ) violate s MIR upper 



m 



Eckart et al. 



limits on Sgr A*, and so does "model ABa<^3' 
(2009 ), in spite of the NIR exponential break invoked to avoid 
the over-production of NIR flux. 

We note also that a similar conclusion, disconnecting the 
NIR and radio portions of the spectrum in the frame of the 
[van der Laan ( [19(56^ model, was reached by Mirabel et al. 
(1998| in the study of discrete ejections from the microquasar 
GRS 1915-H105. 

Expansion models for the sub-mm flares of Sgr A* should 
have two major telltale signatures: spectral indices in the sub- 
mm/mm range indicating an evolution from a thick to thin 
regime over the course of a flare, and high fluxes in the FIR 
(100-300 jjm. Fig. |9] top panel). The upcoming observations 
of the Galactic nucleus with Herschel and the Stratospheric 
Observatory for Infrared Astronomy (SOFIA) will perhaps allow 
the testing of these predictions. 
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Fig. 11. SED of Flare B. L'-, K^-hsnd, and X-ray peak fluxes are plotted in red, orange, and blue, respectively. NIR error bars 
account for the uncertainty in the extinction correction. The 1-100 GeV spectrum of IFGL J 1745. 6-2900 is indicated in dark gray. 
Models SSCl, SSC2, and SB correspond to the green dashed, magenta solid, black dash-dot lines, respectively. 



5.2. SED modeling 



6. Conclusions 



Though sub-mm and NIR flares may be unrelated, NIR and X- 
ray flares clearly derive from the same radiating particles given 
their simultaneity. Past theories of Sgr A* flares usually invoked 
a combination of synchrotron and inverse Compton mechanisms 
(Markoff ltal.pOOlJ |20 07 ; Liu & Melia|2002| |Liu et al.||2004| 
' ' ~ ~ ' 2003 ( |2004| ). We analyze here synchrotron 



2006a b; Yuan et al. 



and SSC emission for Flare B, for which we have NIR spectral 
information. This imposes a relatively hard electron spectrum 
{p < 3, see Sect. |4]l contrary to model I for Flare A, which had 
a softer spectrum (p - 5) as required by the expansion scenario. 
Since we do not have spectral information in X-rays, we pro- 
pose three possible models (see Tab. [5]and Fig. 
the SED of this flare, each having an X-ray si 
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to represent 



lope either hard 
(SSCl), intermediate hai'd (SSC2), or soft (SB). In afl three, the 
NIR flare is due to synchrotron emission. For SSCl and SSC2, 
the X-ray flare results from SSC emission, which is the most 
commonly used configuration in the literature. Nonetheless, if 
the X-ray slope is soft, as for Flare #2 in Porquet et al. (2008 1, 
then the SSC modeling is no longer valid, since it requires an 



enormous magnetic field strength of order 1000 G (Dodds-Eden 
|et al."2009'; 'Trap et al.| [20T0| . In this case, a broadband syn- 
chrotron spectrum with a break (SB) due to radiative cooling 
could give rise to a soft X-ray flare (see Dodds-Eden et al.|2009 
lOTObllTYapet al.|2010[ 



We note in Fig. 1 1 that these flare spectra are all well below 
the 1-100 GeV spectrum of IFGL J1745.6-2900, which is con- 
sistent with the non-detection rep orted in §]3.4| In the pr ogress 
of this work, we noticed a paper by jKusunose & Takahara ( 2010| l 
modeling the time evolution of purely synchrotron (SB) scenar- 
ios for X-ray/NIR flares of Sgr A*. They predicted possible SSC 
counterparts observable with Fermi in a specific class of models, 
where the ratio of magnetic energy over particle energy is small. 
We have shown here that there was no detectable GeV excess 
during Flare B but, we stress that given our low statistics, we 
could only constrain a putative GeV flare to a flux smaller than 
a few times 10"^ erg s"' cm"^, whereas their "model C" induced 
a flux of only ~2 x 10"'' erg s"' cm"^ at maximum. Except for 
an extraordinarily bright event, it is very likely that such models 
will not be testable with Fermi data. 



Only a handful of flaring events have been observed simultane- 
ously accross a broad range of wavelengths up to now. Here, 
we have reported the detection of two new X-ray/NIR flares, a 
strong sub-mm outburst possibly associated with one of them, 
and their GeV y-ray non-detection. These data allowed us to fur- 
ther characterize Sgr A*'s variability and test non-thermal radia- 
tive mechanisms usually employed to describe the physics of the 
flares. 

One flare was followed at the VLT in two NIR spectral bands, 
thus demonstrating the feasibility of a filter switching tech- 
nique to measure the color of Sgr A* with NACO. The result- 
ing spectrum is found to be consistent with a power-law of index 
a - -0.4±0.3 (Fy oc v"), which is a "blue" spectrum (J3 > Q with 
vFy oc v^), implying that NIR flares arise from a particle popula- 
tion different from the one responsible for the quiescent sub-mm 
bump. For this flare, we could also check that there was no de- 
tectable in sync variation of the GeV y-ray GC source. The other 
NIR/X-ray flare was followed by a strong sub-mm flare, and was 
used to test the adiabatic expansion model which is often used 
to account for the longer wavelength flare delays. In this con- 
text, we find it difficult to relate the NIR and sub-mm events by 
the same single synchrotron plasmoid, as it would induce a too 
steep infrared spectrum in regards of other flares. Even assum- 
ing the NIR/X-ray event just acted as a trigger for the plasmon 
expansion, the long delay observed between the NIR/X-ray and 
sub-mm peaks spreads the theoretical light curve of the sub-mm 
flare too much to reproduce the data. However, we suggest an 
expanding plasmon could still drive the sub-mm flare light curve 
provided it started expanding about a couple of hours after the X- 
ray/NIR event. These models can all be tested by radio observa- 
tions obtained in April 2009 and predict a significant flux in the 
FIR that could be observed by the Herschel and SOFIA obser- 
vatories in the near future. Hopefully, the source expansion sce- 
nario will also soon be put to test by VLBI direct size measure- 
ments in the sub-mm range ( Doeleman et al.|2009) . These stud- 
ies of SgrA*'s present variability and the ones p ursuing giant 
outbursts experienced by Sgr A* in a recent past ^Terrier et al. 
2010; Ponti et al. 2010), should thus gradually bring an overall 
characterization of our Galactic nucleus' activity and help us put 
it into the general perspective of AGNs. 
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